Phase Diagram of One-Dimensional Extended Hubbard Model at Half Filling 
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We reexamine the ground-state phase diagram of the one-dimensional half-filled Hubbard model 
with on-site and nearest-neighbor repulsive interactions. We calculate second-order corrections to 
coupling constants in the g-ology to show that the bond-charge-density-wave (BCDW) phase exists 
for weak couplings in between the charge density wave (CDW) and spin density wave (SDW) phases. 
We find that the umklapp scattering of parallel-spin electrons destabilizes the BCDW state and gives 
rise to a bicritical point where the CDW-BCDW and SDW-BCDW continuous-transition lines merge 
into the CDW-SDW first-order transition line. 
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Electronic correlations in solids have been a subject 
of intensive research over the years. Correlation effects 
have the strongest impact at commensurate band filling, 
where a system often undergoes a Mott transition. The 
one-dimensional (ID) extended Hubbard model (EHM) 
with the nearest-neighbor repulsion V, in addition to the 
on-site repulsion U, is a simple, but nontrivial model that 
exhibits rich phase structure The model has a long 
history of research, and considerable amount of knowl- 
edge has been accumulated. Much effort has been de- 
voted to understanding its ground-state phase diagram at 
half filling. In the strong-coupling limit jl], ||, one can 
show that the model has two insulating phases, the spin- 
density-wave (SDW) phase and the charge-density-wave 
(CDW) phase, which are separated by a first-order tran- 
sition line located at U ~ 2V. In the weak-coupling limit 
the perturbative renormalization group (RG) analysis |IJ 
concluded that there is a continuous phase transition be- 
tween the CDW and SDW phases also at U = 2V. It 
was then considered that, as the coupling constants in- 
crease, the continuous-transition line changes into the 
first-order one at a tricritical point in the intermediate- 
coupling regime. This picture was supported by both 
numerical||, [|, || and analytical H 0] studies and had 
been regarded as the complete phase diagram of the EHM 
at half filling. 

Quite recently, however, NakamuraQ found numeri- 
cally that another phase exists between the CDW and 
SDW phases for weak couplings. The new phase is the 
bond-charge-density-wave (BCDW) phase in which the 
Peierls dimerization occurs spontaneously. He concluded 
that SDW-BCDW and BCDW-CDW transitions are con- 
tinuous and that these two transition lines merge at a 
multicritical point into the first-order line separating the 
CDW and SDW phases 0. His claim was confirmed by 
a recent extensive Monte Carlo calculation pi. The ap- 
pearance of a spontaneously dimerized phase in the EHM 
is surprising and calls for thorough theoretical study. So 
far the BCDW phase has been shown analytically to ex- 
ist only in models with extra correlated-hopping interac- 
tions For the original EHM, however, the origin of 



the BCDW phase and the nature of the associated phase 
transitions are not fully understood. In this Letter, we 
will provide theoretical argument for the existence of the 
BCDW phase by reformulating the weak-coupling theory 
to include higher-order terms. Using the bosonization 
technique, we derive a set of RG equations and discuss 
the critical properties of the phase transitions. We find 
that the umklapp scattering between electrons with par- 
allel spins is responsible for the emergence of the bicritical 
point. 

The Hamiltonian of the ID EHM is 
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where n^ a = ct ff Cj )<T 



and 



denotes the creation operator of an electron at the jth 
site with spin a. Following the previous studies on mod- 
els with correlated-hopping interactions ]ll| , we consider 
the CDW, SDW, BCDW and bond-spin-density-wave 
(BSDW) phases. They are characterized by the order pa- 
rameters, Ocdw = (— l) J 'rij, Osdw = (— l) J '( n i,T ~ n i-\)-> 



Obcdw = Ea(' 
(-l) J '(ct t C, ul t - ct , c„- 
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Wc first focus on the weak-coupling limit U, V -C t. 
The hopping t generates the energy band with disper- 
sion £k — — 2icosfc. At half filling the Fermi points 
are at k = ±kp = ±7r/2a, where a is a lattice con- 
stant. Electrons experience two-particle scattering by 
the on-site and nearest-neighbor repulsions U and V. 
We follow the standard g-ology approach ]]], [l3| and 
parametrize the scattering matrix elements by the cou- 
pling constants g. In lowest order in U and V they are 
known to be = —2Va, gi± = (U — 2V)a, ffail = 2Va, 
92X = (U + 2V)a, 5 3|| = -2Va, g 3 ±_ = (U - 2V)a, 
<74|i = 2Va, Qix = (U + 2V)a, where we have used 
the standard notation [13) . Here we note that both 
gi± and g%± vanish at U = 2V, and this is the rea- 
son why the lowest-order calculation predicts the direct 




FIG. 1: Vertex diagrams. Solid lines denote the low-energy 
states and the dashed lines represent the high-energy states 
to be integrated out. 



CDW-SDW transition at U = 2V. Hence, we need to 
go beyond the lowest order to see if the BCDW phase 
really exists. To this end, we adapt the two-step RG 
scheme used in Ref. Q . (i) We separate the states into 
low-energy states (||fc| — kp\ < A) and high-energy ones 
(||fc| — &f| > A) by introducing a momentum cutoff A, 
and integrate out high-energy states to obtain effective 
scattering matrix elements for low-energy states, (ii) We 
then derive one-loop RG equations for these matrix el- 
ements using the standard bosonization method. The 
diagrams for the effective couplings up to second order 
in U and V are shown in Fig. |[ The explicit calculation 
yields 



9i± = (U-2V)a 
g 3± = (U-2V)a 
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where Ci(A) = 21n[cot(aA/2)] > and C 2 (A) = 
2cos(aA) > 0. The weak dependences of C^s on A al- 
low us to set A = 7r/4a; different choices will only lead 
to small quantitative changes. We see that gi± < and 
.93± > at U — 2V due to the Ci term. This implies 
that a new phase different from the CDW and SDW can 
appear for U ~ 2V, as we will show shortly. The zeros of 
gi± and g 3 ± are shifted from U = 2V due to the momen- 
tum dependence of the matrix element 2Va cos(qa) for 
the virtual scattering of high-energy states (q ^ 0, 2kp). 
There is no symmetry principle that enforces gi± and 
g 3 ± to vanish simultaneously. 

Having derived the effective scattering matrix elements 
for low-energy states, we now apply the bosonization 
method. The right-going and left-going electron fields 
^±,<t are written jj], [b| 



: exp [ipk F x + ip tp p ^(x)] , (4) 



where (p p . a (p — +/ — ) are the chiral bosonic fields and 
\j]u , T]cr' } = 2(5 CTi(T /. The bosonic fields obey the commu- 
tation relations [ip p , a (x), <p P ,a' {x')] = ipirsgn(x — x') 8 a ^i 
and [ip+ t0 -(x), <p- l(T > (x 1 )] = iirSa-.o-i. We define chiral 
charge fields, 9 p = (<p p ,t + tp p ^)/2, and chiral spin fields, 
4>p = {fpA ~ Vp,l)/2, to write the Hamiltonian density 
for low-energy states: 



H = 77- Y] [v P (d x 9 p f + v a (d x (j>p) 
Ztt z — ' 
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FIG. 2: Phase diagram obtained by minimizing the potential 
energy. Bicritical points are at g c ~ —gs = ±<?cs. 
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a 2 {dj+){d x 9-){d x (t> + ){d x (j>-), (5) 



The renormal- 
94± - 3i||)/27r and 



p=+,- 



where 9 = 9 + + 6>_, 4> = 4> + 4 
ized velocities are v p = 2ta + (<? 4 | 
v a = 2ta+((74|| — <?4_l — 5i || ) /27r. To simplify the notation, 
we have written g c = g 3± , g s = gu_, and g cs = g 3 \\. The 
other coupling constants are given by g p = g2_L + <?2|| ~9i\\ > 
9a = 92± - 92\\ + 9i\\, and g ps = g ca = g pa = -2Va to 
lowest order in V. The g ps (g pa ) coupling comes from 
the backward scattering of electrons with opposite (par- 
allel) spins, while the g ca coupling is generated from the 
umklapp scattering of electrons with antiparallel spins. 
The SU(2) symmetry in the spin sector ensures g a = g s , 
gcs = gca, and g ps — g pa , and therefore it is important 
to retain the g pa term. 

In terms of the phase fields 9 and <fi the order param- 
eters are written as 

Osnw(x) °f cos9(x) smtfi(x), (6) 

CcdwC^) °c sin9(x) cos(f>(x), (7) 

Cbcdw^) oc cos9(x) cos0(x), (8) 

C , BSDw(a ; ) oc sin6*(a;) sm<p(x). (9) 

The phase diagram can be qualitatively understood via 
a quasi-classical analysis: we neglect spatial variations 
of the fields and focus on the potential, V(9, </>) = 
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— g c cos 29 + g s cos 2(f) — g cs cos 29 cos 2(f>, where g cs = 
53 1| < 0. The order parameters take maximum ampli- 
tudes when the fields 9 and <fi are pinned at the follow- 
ing potential minima: (9,4>) = (0, ±7r/2) in the SDW 
state, (±7r/2,0) in the CDW state, (0,0) or (tt,0) in 
the BCDW state, and (n/2,±ir/2) in the BSDW state 
(mod 7r). In these states the potential energy V(9, <j>) be- 
comes VsDW = ~9c ~ 9 s - \g cs \, V" C DW = 9c + 9s - Iffcsl, 

Vbcdw = ~g c + g s + \g cs \ and Vbsdw = g c - g s + \g cs \, 

respectively. Comparing these energies, we obtain the 
phase diagram in the g c -gs plane (Fig. 2). The direct 
CDW-SDW transition is first order because there is a 
potential barrier of height min(|gr cs |, 2|g cs | — 2|g c |) be- 
tween the corresponding minima. The other boundaries 
located at g s = ±|<? C s| an< 4 9c = ±|<7cs| are continuous 
transitions, because the pinning potential for 9 or <f> van- 
ishes when the other phase field is pinned. When g cs = 0, 
the first-order CDW-SDW transition line collapses to a 
tetracritical point. 

To obtain the ground-state phase diagram of the EHM, 
we need to include the renormalization of the coupling 
constants due to quantum fluctuations of the fields. A 
systematic analysis in the weak-coupling limit can be 
done by applying the perturbative RG method to H 
(||). The SU(2) spin symmetry guarantees the relations 
g a = g s , g ca = g cs , and g pa = g ps to hold in the scal- 
ing procedure. The one-loop RG equations that describe 
changes of the coupling constants during the scaling of 
the short-distance cutoff (a — > ae dl ) are then given by 



d 

dl Gp - 


+ 2 G 2 + G 2 S + G s G ps , 


(10) 


d G 

dl Gc ~ 


+ 2 Gp G c — G s G cs — G cs Gp S , 


(11) 


d G - 
dl Gs ~ 




(12) 


d G 


— 2 G cs + 2 G p G cs — 4 G s G cs 





— 2 G c G s 2 G c Gps — 4 G cs G ps , (13) 



^■Gp S — — 2 Gps + 2 G p G s 

-4G c G cs -4G 2 cs -4G s G ps , (14) 

where G v = g v /{Awta). From Eqs. (|||) and @ one 
finds that g cs and g ps are irrelevant and renormalized to- 
wards zero for weak interactions. It is therefore natural 
to ignore g cs and g ps first. With this approximation the 
Hamiltonian reduces to two decoupled sine-Gordon mod- 
els, and it is easy to follow the RG flows of G p , G c , and 
G s from Eqs. Since g p = (U + 6V> > 0, G c 

is relevant and grows at low energies. The coupling G s is 
marginally relevant (irrelevant) for g s < (<? s > 0). The 
phase diagram of the EHM is obtained tentatively from 
Fig. 2 by setting g c = g 3 ± [Eq. ©], g s = 9xx [Eq. @], 
and g cs = 0. When U is sufficiently larger than 2V such 
that g c > and g s > 0, we have the SDW phase. If U is 



smaller than 2V (g c < and g s < 0), then we have the 
CDW phase. Around the U = 2V line we find the BCDW 
phase, where g±± < and g^i_ > due to the Gi term. 
The BSDW phase does not exist in the EHM. The charge 
excitations are gapful except on the CDW-BCDW tran- 
sition line where the relevant pinning potential vanishes. 
The spin excitations are gapless in the SDW phase and 
on the SDW-BCDW transition line and gapful otherwise. 
In the gapped phases the charge gap A c and the spin gap 
A s are given by A c ~ t\G c \ 1/2G » and A s ~ £exp(l/2G s ) 
for |G C | <§; 1 and Q < -G s < 1, respectively. 

Next we examine effects of the parallel-spin umklapp 
scattering g cs for U ~ 2V. Let us assume U — 2V — 
-C 2 V 2 /irt + Oiy z lt 2 ), i.e., g c w and g s < 0. We are 
considering the situation very close to the CDW-BCDW 
transition. In this case the spin gap is formed first as the 
energy scale is lowered, and we can replace cos2<^ with 
its average (cos2</>) ~ (A s /t) 2 for energies below the spin 
gap. This means that the cos 29 potential that tries to 
pin the fluctuating 9 field has the effective coupling 

9* c = 9c + g cs (cos 2(f>). (15) 

The CDW-BCDW transition occurs when g* = 0, i.e., 
g c = —gcs (cos 2cf>) > 0. The phase space of the 
BCDW state is reduced upon inclusion of the g cs term. 
Note, however, that the CDW-BCDW boundary does not 
move across the U = 2V line because |g cs (cos 2<f>) \ ~ 
2Va exp[— c(t/V) 2 ] is much smaller than the G2 term for 
V <C t, where c is a positive constant. A similar argument 
applies to the region near the SDW-BCDW transition. 
Suppose that U - 2V = +C 2 V 2 /irt + 0{V 3 /t 2 ) (g s « 
and g c > 0). In this C&1S6, (IS the energy scale is lowered, 
the charge gap opens first and the 9 field is pinned at 
9 = (mod 7r) . Below the charge-gap energy scale the </> 
field is subject to the pinning potential g* cos 2<j> with 

9* s = 9s -g cs (cos29), (16) 

where (cos 2(9) ~ (A c /t) 2( - 1 ~ G »'> . The SDW-BCDW tran- 
sition now happens at g s = g cs (cos 2(f)) < 0, and thus the 
SDW-BCDW transition line moves to increase the SDW 
phase. Again the phase boundary is not changed beyond 
the U = 2V line as |# cs (cos20)| ~ 2Va{c'V/t) 7Tt ^ v is 
much smaller than V 2 a/t, where d is a constant of or- 
der 1. This completes our proof of the existence of the 
BCDW phase near the U = 2V line in the weak-coupling 
limit. 

For larger U and V, the g cs coupling becomes less irrel- 
evant, and the BCDW phase will eventually disappear. 
Since the cosine factor in Eqs. @ and @ can be con- 
sidered as renormalization of g cs , we conclude that the 
two continuous lines meet when the renormalized cou- 
plings satisfy the relation 

G c = -G s = -Gcs = G (G > 0) (17) 



4 



4 



ft 



CDW 




BCDW 


SDW 







4 



U/t 



8 



FIG. 3: Phase digram of the half-filled extended Hubbard 
model. The bicritical point is at (U C ,V C ) — (5.0t,2.3t). 

in the low-energy limit. This is the condition for the 
bicritical point in the RG scheme. Note that the con- 
dition is not simply that the g cs term becomes relevant, 
as previously assumed ||. When Eq. ( |l7| ) is satisfied, 
the effective potential takes a simple form V(9, 4>) = 
—G(cos29 + cos 20 — cos 29 cos 20), which has an inter- 
esting feature that its potential minima are not isolated 
points but the crossing lines 9 — nm or <p = 7rn (m, n: 
integer). On these lines either 9 or <j> becomes a free field; 
the theory has more freedom than a single free bosonic 
field, but less than two free bosonic fields. We thus ex- 
pect that the theory of the bicritical point should have a 
central charge larger than 1 but smaller than 2. Detailed 
analysis of the critical theory is left for a future study. 

We have numerically solved the scaling equations (|l0|)- 
( |l4| ) to obtain the global phase diagram of the EHM. The 
phase is determined by looking at which of the couplings 
G c , G s , and G cs becomes relevant. The idea is essen- 
tially the same as what we have discussed above. If \G C \ 
grows with increasing I and reaches, say, 1 first among 
the 3 couplings, then we stop the integration and cal- 
culate G* = G s — G cs . Since the charge fluctuations 
are suppressed below this energy scale, we are left with 
Eq. (|T|) where G s replaced by G* and G cs = 0. We im- 
mediately see that a positive (negative) G* leads to the 
SDW (BCDW) state. If \G S \ becomes I first, then the 
sign of G* = G c — G cs determines the phase: the CDW 
(BCDW) state for G* < (G* > 0). Finally, when 
| G c ,j | reaches 1 first, we stop the calculation and com- 
pare 



marginally irrelevant, and therefore the spin sector be- 
comes gapless. The phase diagram obtained in this way 
is shown in Fig. ||. For weak couplings the BCDW phase 
appears at U ~ 2V , and the successive continuous tran- 
sitions between the SDW, BCDW and CDW states oc- 
cur as V/U increases. As U and V increase along the 
line U — 2V, the BCDW phase first expands and then 
shrinks up to the bicritical point (U c , V c ) w (5.0t, 2.3t) 
where the two continuous-transition lines meet. Beyond 
this point the BCDW phase disappears and we have the 
direct first-order transition between the CDW and the 
SDW phases. The phase diagram (Fig. ^) is similar to 
the one reported recently || E3. The position of the 
phase boundaries in Fig. || is not reliable quantitatively 
however, as we have used the perturbative RG equations 
which are valid only in the weak-coupling regime. 

In summary, we have studied the ground-state phase 
diagram of the ID extended Hubbard model with repul- 
sive interaction at half filling. We have shown analyti- 
cally that the BCDW phase appears at U ~ 2V in the 
weak- interaction limit. We have also discussed the in- 
stability of the BCDW state and the emergence of the 
bicritical point due to the parallel-spin umklapp scatter- 
ing. 
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are already suppressed by the G cs eus20eus20 puleiilial, 

we can deduce the phase from the quasi-classical argu- 
ment. From Fig. 2 we see that we have the SDW state 
for G s > -G c and the CDW state for G s < -G c . In 
the SDW state the pinning potential for the <j> field is 
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